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Abstract 

Analytical expressions are given for the eigenvalues and eigenvectors of a Hamilto- 
nian with su g (2) dynamical symmetry. The relevance of such an operator in Quantum 
Optics is discussed. As an application, the ground state energy in the Dicke model is 
studied through su 9 (2) perturbation theory. 
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1 Introduction 



Many models in Quantum Optics, such as Raman and Brillouin scattering, parametric 
conversion and the interaction of two-level atoms with a single-mode radiation field (Dicke 
model), can be described by interaction Hamiltonians of the form (see, e.g. [1-4]) 



H 



( Ai ... \ 

A x M-\ ••• 

... A^i +2 A_i +1 

V ... A_ l+1 J 



(1.1) 



The dimension d = (21 + 1) of this matrix is by no means small (for instance, in the Dicke 
model, 21 is just the number of atoms considered). Therefore, the finding of analytical 
expressions for the corresponding eigenvalues and eigenvectors of H is essential in order 
to solve the dynamics of the model. It is also important to point out that, in some cases, 
H can be seen as a perturbation of the J x generator of an underlying su(2) dynamical 
symmetry. This fact has been successfully used in order to describe many features of these 
models [3], and it will be also relevant in what follows. 

In this paper we show that the su 9 (2) quantum algebra (see e.g. [5-13]) can be used to 
define a Hamiltonian of the type (1.1) in a natural way. Such Hamiltonian is introduced as 
a simple function of the su g (2) generators having as non-deformed limit the J x generator 
of the su(2) algebra. By considering the well-known representation theory of su g (2) (which 
is revisited in Section 2), the Hamiltonian is defined and its eigenvalues and eigenvectors 
are found (Section 3). The spectrum obtained is essentially anharmonic; thus we have a 
new exactly solvable nonlinear quantum model with su g (2) dynamical symmetry. 

In Section 4, relevant Clebsch-Gordan coefficients for both the bare and the dressed 
basis are considered. Finally, in Section 5 we present an application of the previous 
analytical results to finding the energy of the ground state of the Dicke model by making 
use of a perturbation theory around the su g (2) Hamiltonian. This preliminary study shows 
that the (/-deformed J x operator here introduced is physically meaningful in the context of 
the quantum optical Hamiltonians mentioned above (compare with the results given in e.g. 
[2-4]), and the explicit dynamical features of the su g (2) Dicke model will be fully developed 
in a forthcoming paper. In this work, we shall provide the basic algebraic properties that 
are needed in order to solve it explicitly. It is interesting to stress that some of these 
properties show new (to our knowledge) features of the su g (2) algebra, all of them related 
to the deformed J x = (J+ + J-)/2 generator. 

2 The su^) algebra 

Let the operators J z , J± generate the quantum algebra su g (2) with the coproduct [5]-[13] 

A(4) = 4®g" J Hg Jz ®J±, (2.1) 
A(J*) = J z ® 1 + 1 ® J z . (2.2) 
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Deformed commutation rules consistent with the previous map are given in the form: 

[J Z ,J±] = ±J±, [J+,J-] = [2J Z ], (2.3) 

where [x] := q ~ q _i and q = e z l 2 . We shall assume that q is not a root of unity, and we 
shall recover "classical" results when q — > 1. 

Let us introduce, firstly, the "bare" basis of eigenvectors of J z , 

2J z \l,m) = 2m\l,m). (2.4) 

In this basis, the (2 I + l)-dimensional irreducible representation of su g (2) is given by (2.4) 
and 

J±\l,m) = y/[lTm][l±m + l]\l,m + l). (2.5) 
Throughout the paper we will use the following (standard) notation for (/-numbers, 

[n] = q~ n+1 + q~ n+3 + q- n+5 . . . + q n ~ x = q _ J_ 1 , [0] = 0, [-n] = - [n] . (2.6) 
We shall also introduce the symbols 

[n, m] := q'^ m + q '^) m + q ~(n-5)m _ _ _ + g (n-l)m = [ nm ]/[ m ] . ( 2 .7) 

Thus, lim g _^i[n, k] = n, and 

[n,l] = [n], [fc,0]=A;, [2, fc] = q k + q~ k = [2k] /[k] . (2.8) 

3 The Hamiltonian 

Let us consider the following su g (2) operator as a Hamiltonian: 

H = q J */ 2 (J + + J_)q J >/ 2 . (3.1) 

We stress that the q — > 1 limit of H is just 2 J x , but is not the standard way to introduce 
the deformed J x operator in su g (2). On the other hand, one can easily check that, from 
the definition (3.1) and the maps (2.1)-(2.2), the coproduct of H can be deduced (note 
that A is an algebra homomorphism: A(XY) = A(X) A(Y)): 

A(H) = H®l + q 2J *®H. (3.2) 

Such a form of the coproduct is known (see, e.g., [13]). In the (21 + l)-dimensional 
representation Di, H thus takes the form (1.1) with 

A m = q m - 1/2 ^[l+m][l-m + l}. (3.3) 
In particular, when / = 1 we have, 

/ q 1/2 VW\ \ 

D^H) = q^ 2 V\2] q~ 1/2 VW] • (3-4) 

V o q-^vm J 
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A straightforward computation shows that the spectrum of this operator is [2], 0, —[2]. 
The corresponding normalized eigenvectors are 



1 



,1/2 



a- 1 / 2 





-Q i/2 vm I 



(3.5) 



Due to the richness of the structures underlying quantum deformations, this remarkable 
anharmonic deformation of the J x Hamiltonian can be explicitly solved for arbitrary / as 
follows. 



3.1 Spectrum and eigenvectors 

It can be proven (by induction and using the coproduct to construct higher dimensional 
representations from the three-dimensional one considered before) that the spectrum of 
this operator for a given / is just [2m], with m = —I, . . . , I. 

Moreover, the eigenvector \l,l)' corresponding to the highest eigenvalue [2/] is given 
by the following components in the bare basis \l,m): 

a ml ee (l,m\l,l)' = 9 -( J -V2) LZ B_Z. (3.6) 
— y [t+mj![/ — my. 

Here, the prime means that a special normalization is accepted, temporarily, where an = 
ql(l-i/2) _ i^ a _ [ _ lm xhe next eigenvector \l, I — 1)' with the eigenvalue [21 — 2] is given by 



cw-i = (l,m\ l,l- 1 )' = Qhm q- 2m+ \fW] (l - q 1 ^' 1 [2] Lj^) . (3.7) 
Finally, an arbitrary eigenvector with eigenvalue [2m] can be deduced, namely 



a mn = (I, m|U)' = a ml ^ _JL (3 . 8 ) 

x V7-1 y„-ia-™-2n)+i(i+i)/2 [2(Z-n)]!! [/-m][/-m-l]...[/-m-j + l] 

^ [j]![2(/-n-j)]!! [2Z][2Z-1]... [2Z-J + 1] ' { -> 

where [2n]!! := [2n] ■ [2n - 2] ■ ■ ■ [2]. 
3.2 Normalization 

These eigenvectors can be easily normalized in terms of (/-numbers. For the ground and 
top states we have, 

|M) = |M)' , where Aft = A^ = '(IJjlJ)' = q' 1 ^ g fj^pZ^T ( 3 - 10 ) 
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A simple calculation gives, 

Kl = I[V + ?"*) = 2 K (l 2 + T 2 ) • • • • • {q 2l - l +q- 2l+l ) = n' [2, k] . (3.11) 

fc=0 fc=0 

The excited states can be normalized as follows: 

l+n -i l+n 



" " k=—l+n ' J k=—l+n 

For instance, 

A^ 1 2 /2 , 1/2 =AA 1 2 A _ 1/2 = 2, AT 1 2 , 1 =AT 1 2 ; _ 1 = 2[2], AA 2 = [2] 2 , ... (3.13) 

Note that all of the coefficients A/" 2 n go to 2 2 * in the limit q — > 1. Finally, if we denote the 
entries of the normalized eigenvectors as 

A l mn = a l mn /M ln , (3.14) 

the following relations hold, 

A l mn (q) = (-l) l - n A l _ m>n (l/Q) , A l mn (q) = (-iy- m <_„(<7) . (3.15) 

They generalize the known symmetry of the non-deformed su(2) case (q = 1) by involving 
the transformation q q^ 1 . 

3.3 Orthogonality relations 

For physical applications (like the calculation of mean values) is desirable to further develop 
the previous "(/-arithmetics" . With this goal in mind, let us note that the quantities \cx m i\ 
from Eq. (3.6) play the role of the binomial coefficients. In particular, the moments of 
this deformed binomial distribution, i.e. the mean values, 

21 

<Y k > 2l = Y J Wk-iA 2 Y k . (3.16) 

k=0 

can be calculated. With this definition, the equation (3.10) is then rewritten in the form, 

21-1 

<i>2i= n(9 s +<r s ). (3.17) 

s=0 

Moreover, one has the following deformed formulas (0 < j,M < 21) for the binomial 
distribution: 

2l-j-l 

<q~ 2jk >2l= q~ 2jl IT + (3-18) 
s=-j 

< [k][k-l] . . . [ k -M + l]q- k( V +M "> > 2l = q -2Mj-M{M+l)/2 < g -2jk >a _ M (319) 

2l-M-l-j 

= [2l)[2l~l)...[2l-M + l]q~^ 2l+M ^- M ^ M+1 ^ 2 H (q s + q' 3 ) . (3.20) 

s=-j 

These formulas are equivalent to the fact that the dressed vectors form an orthonormal 
basis. 
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4 Clebsch-Gordan coefficients for the dressed basis 



We can decompose the tensor product of the dressed vectors (3. 9), (3. 14) into the irre- 
ducible parts: 

j=l+i 

\Lm) ® \h±) = E Cli,mi\j' m + *)i2 , (4-1) 

3=1-1 

where Cj lmi are the Clebsch-Gordan (C.G.) coefficients for the dressed basis. However, 
the vector \l, m) ®\l,i) in the left-hand side and vectors \j, m + i)\2 in the right-hand side 
are eigenvectors of J_ x in the corresponding representations only for i = 0. Thus, only the 
coefficients mQ are of interest. For 1 = 1, they are given as 



£11,10 = y j^/j ' £11,10 = -\J q i +q --i > £11,10 = °> I 4 - 2 ) 
£11,00 = \/ foi roi ' £11,00 = 0' £11,00 = — i/t^p (4-3) 



[2] [3]' 



£11,-10 — £11,10 > £11,-10 — ~~ £11, ioj £11,-10 — 0- (4-4) 
They thus differ from the coefficients in the bare basis [9] (though have the same limit 

By rewriting the definition, 

£/l;p,0= 12<J 1 p|i 1 P>l|Vq>2, (4.5) 

in terms of components of the dressed vectors, and afterwards replacing explicitly the 
components of 1 1,0)2, we get the relation 

£^o = E < (^ 1/2 ^-1,1 4-i,p - ^+1,-1 4 + i,p) • (4-6) 



fc=-jf 

Here, ^.^n = if |m| > L and 

C'^mi.ma = Mj,™! + m 2 \l,m 1 ) 1 \l,m 2 )2 , (4.7) 

are C.G. coefficients in the bare basis. Another useful formula for the Clebsch-Gordan 
coefficients in the dressed basis can be obtained as follows. Starting with the expansion 
of the tensor product, 

|/^) 1 |l 1 0) 2 = E£n ;P oliP)i2, (4-8) 

3 

we can rewrite the dressed vectors in the bare basis, \j,p)i2 = J2in=-j •A- > m , p \j-, m)i2, and 
the vectors \j, m)\ 2 in terms of the tensor product, using bare C.G. coefficients, 

1 

\j,m) 12 = C/i ;m - m2 , m2 l^™-™2>i|l,™ 2 }2. (4.9) 

m,2=— 1 
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We have, 

A mip A m2 Q = X] —ll;p0^m 1 +m 2 ,p^ll:m 1 m 2 - (4-10) 
3=1-1 

Multiplying by Cf i , making a summation under the condition m,\ +m 2 = m, using 

the orthogonality of the Clebsch-Gordan coefficients, 

X] < ^hl2-,m-m 2 ,m 2 ( ^hl2\m-m,2,m2 = ^3 k ' (4-H) 

and, finally, replacing the explicit form of A^q, we have the connection between the 
coefficients in the bare and dressed basis, 

C j A j - 1 f A 1 -, o 1 / 2 ^' + A 1 , cr^ 2 C j \ (4 12) 

—j<m<j, — I < p < I . 

For instance, using these formulas for m = j, we can express the "dressed" coefficients in 
terms of "bare" ones. 



Finding components of the first three rows of the matrix A mn , 



A 



r>i i \ \ \ kk i.K >\\ r, '■■ > i 1..1.JA: nidi i i a w ' / 
n 2 -//2 



[2Z]! 



V [Z-n]![Z + n]! ' 



g n»-(3J-l)/2 r [2Z ], - [ 2 „] 



" M» V P — r*]![«H-r.]! v/pj' 1 j 

, g" 2 -^ 2+1 / [2/]! ~ q 1 - 2l [2n] 2 -[21] 

l ' 2 ' n ~ Mn V [Z-n]![/ + n]! J[2l] [21] [21 - 1] ' 

and, from (4.12) we arrive at the explicit expressions, 



'[Z + l+p][Z + l-p][2,Z + l+p][2,/ + l-p] 



^l;p0 



[2][2Z + 2][2Z + 1] 

[2p] 

V[2Z + 2j[2Zj ' 
'[Z+p][Z-p][2,Z+p][2,Z-p] 



[2][2Z][2Z + 1] 



5 Application to the Dicke model 



In order to demonstrate how our approach works we will apply it to the Dicke model, 
which describes the interaction of a system oi N = 21 two-level atoms with the quantum 
radiation field in an ideal cavity. (This model is mathematically equivalent to the three- 
photon Hamiltonian describing three-wave mixing.) The Hamiltonian can be written in 
the matrix form (1.1) with the matrix elements 

B m = yj(l + m)(l + 1 - m){s + 1 - m), (5.1) 
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where s > I is a parameter (2s is the excitation number which is conserved in this model 
[3]). We restrict ourselves with the case of the highest nonlinearity, s = I. (In the language 
of the three- wave mixing processes this corresponds to the second harmonics generation.) 
We are specially interested in the limit of large /, which corresponds to high photon 
numbers. 



We will take the Hamiltonian 



(5.2) 



with H (3.1) belonging to su q {2) as a zeroth-order Hamiltonian and we will find the Dicke 
spectrum by using the perturbation theory. In the present work we restrict ourselves with 
the energy of the ground state, for the sake of simplicity. The values of q and $7 are to 
be chosen. The simplest way to fit them is to provide the coincidence of the points where 
the matrix elements B m of the three- wave Hamiltonian and the matrix elements 



A m = ft q m ~ 1/2 ^[l -m][l + m+l] = n(l,m+ 1\2 J x \l, m) 
of the Hamiltonian Hq take their maximum values. It gives (for s = I) 

1 



3 n v 7 ^- 
a = Nlogq = - log — — 



-0.7218, 



(5.3) 



(5.4) 



and the maxima of B(m) and A(m) occur in the point mo = — (/ — l)/3. We chose the 
coefficient fi to make equal the values of A m and B m in their maxima. This gives 



n 



4(AT+ 1)3/2 



27[N + 1] 



(5.5) 



At the next step, we find the approximation for the three- wave Hamiltonian in the 
form: 



B n 



QA m <p(m), 4>(m) = 1 + faA - 2 A 2 + 3 A 3 



m — mo . 



(5.6) 



We thus restrict the expansion up to the third-order polynomial <p(m). We can find 
explicitly the coefficients 4>j by comparing the Taylor expansions for the matrix elements 
B(m) and A(m) around the point m . We have: 



B(m) 



2(N + If' 2 



and 



A(m) = 



[N + l] 



27 



1 



1 - 



27 



A 



2a 2 



8 VA^ + l 



A 



2 27 
+ 



A 



8 VN + l 



+ 



tanh^a V^ + l 



4a 3 



tanh 2 a V^ + l 



A 



O (iV" 4 ) 



+ N 



which determines the polynomial </>(A): 
2a 2 \ / A ^ 2 



,22 



tanh^a/ V^+l 



27 
+ hr + 



4a 3 



(5.7) 



(5.8) 



8 tanh a 
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Now we may substitute A = m — mo = J z + (I — l)/3 and rewrite (5.6) in the matrix form: 



H » fi [J+cp(J z - m ) + <j>(J g - m )J-} = 2tl {J x , f(J z )} . (5.10) 

Here J± jZ are generators of su q {2) and {A, B} = AB + BA. The new function f(J z ) is 
also a polynomial of degree three, whose coefficients can be easily found. Now the ground 
state energy is approximately given as 

3 

(-i,i\H\-i,i) « -fipz] £ fk(dA\( J z) k \zhl)- ( 5 - n ) 

fc=0 

Therefore, we have reduced the problem to the calculation of averages of the powers of 
the operators J z (the moments) in the eigenstates of the operator J x . 

Though this problem can be solved for arbitrary eigenstates, here we consider the 
ground states only. By using the results of the previous sections we can write the generating 
function for these moments. For arbitrary p = exp [i and q = expi we find: 

Now, let us introduce the notations 

N-l .jv_i 



Sk 



f"- 1 1 r n 

= J2 tanhk J t - dxtanh k xt+ - \tanh k t + tanh k (N - l)tj , (5.13) 

i=i Jl 2 



where we have used the Euler-Maclaurin summation formula. Differentiating (5.12) with 
respect to \i we find, 

(-l,l\2J z \-l,l) = S u (5.14) 

(-M|(2J 2 ) 2 |-M) = Sj-S 2 + N, (5.15) 

{-l,l\(2J z f\-l,l) = Sf-3S 2 S 1 + (3N-2)S 1 + 2S 3 . (5.16) 

On the other hand, for the first three sums we have from (5.13), 

N , , . , tanh(a — t) a , 

51 » — log cosh(a - t + ^ ^ , t = — = log q, (5.17) 

a 2 TV 

5 2 « jV _ 1 _ jV tanh(a-t) + tanh 2 (q-t) ) 

a 2 

TV, w N ,tanh 2 (a-t) tanh 3 (a-t) 

53 w — log cosh(a — t) — N + 1 . (5.19) 



The combination of all these formulas gives the approximation for the energy of the 
ground state. Comparing with the numerical results, we can say that the accuracy for the 
energy of the ground state is 1.5% for 100 atoms (N = 100) and 0.35% for 400 atoms. 
Note that it is the Maclaurin summation formula (5.13) that reduces the accuracy, which 
would be otherwise much higher. However, it gives the correct asymptotic behaviour when 
N — > 00, which is sufficient for our goals. We may mention also that our method produces 
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much better accuracy than the analogous perturbation theory with common su(2) as a 
dynamical symmetry algebra [3, 14] or than the variational method with the su(2) coherent 
states as probe states [15]. The complete description of the spectrum and the dynamical 
analysis of the Dicke model by means of the present approach will be given elsewhere. 
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